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ABSTRACT 

In this paper, we identify Cen A and M87, two nearby FanarofF-Riley I (FRI) radio galaxies, with 
high-energy-peaked BL Lac-hke (HBL-like) objects by investigating their spectral energy distributions 
(SEDs). The SED peak of Cen A at ~150 keV, which was generally believed to be the peak of inverse- 
Compton emission as in the case of 3C 273, is found to be actually the peak of synchrotron emission. 
The synchrotron emission of M87 peaks in the far UV band. We summarize the properties of 7-ray-loud 
blazars, especially those of TeV BL Lac objects, and generalize them to HBL-like FRI radio galaxies 
according to the unified scheme of BL Lac objects and FRI radio galaxies. We infer that Cen A may 
have a peak in its Compton component power output at ~1 TeV, and that M87 may have a Compton 
emission peak at ~0.1 TeV. For Cen A, the estimated TeV 7-ray flux during outburst is F(0.25 — 30 
TeV) = 6.4x10"^ erg cm'^s"!, and for M87, F(0.25 - 30 TeV) = 1.1x10"" erg cm'^s"!. Both 
fluxes are detectable by TeV detectors available today, and hence Cen A and M87 are TeV 7-ray source 
candidates. By investigating the long term variability, we predict that Cen A will undergo an outburst 
in the near future and will be detectable at the TeV 7-ray energy range using the CANGAROO and the 
German-French-Italian experiment HESS TeV 7-ray telescopes. 



Subject headings: BL Lac objects: general 
mechanism: nonthermal - 

1. INTRODUCTION 



- galaxies: individual (Cen A and M87) - 
X-rays: galaxies — gamma rays: theory 
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Up to now, five active galactic nuclei (AGNs) in to- 
tal have been discovered to be TeV {E >0.3 TeV) 7-ray 
sources. All of them are nearby BL Lac objects, namely 
Mrk 421 (z=0.031. Punch et al. 1992), Mrk 501 (z=0.034, 
Quinn et al. 1996), lES 2344-1-514 (z=0.044, Catanese 
et al. 1998), PKS 2155 - 304 (z=0.117, Chadwick et al. 
1999), and 3C 66A (z=0.444, Nespher et al. 1998). BL Lac 
objects and flat spectrum radio quasars (FSRQs) consti- 
tute a rare and extreme blazar class of AGN which is char- 
acterized by strong and rapid variability, high polarization, 
and apparent superluminal motion. These extreme prop- 
erties are generally interpreted as a consequence of non- 
thermal emission from a relativistic jet oriented close to 
the line of sight (Blandford & Konigl 1979). According to 
unified schemes, BL Lac objects are intrinsically the same 
as Fanaroff- Riley class I (FRI, Fanaroff & Riley 1974) radio 
galaxies, while FSRQs are the same as FRII radio galaxies 
(e.g. Urry & Padovani 1995). 

Besides the above five TeV 7-ray-loud blazars, 66 blazars 
have been detected to date as GeV 7-ray sources by the 
EGRET experiment on board CGRO (Mukherjee et al. 
1997; Mattox et al. 1997; Hartman et al. 1999; Sreeku- 
mar et al. 1999). The 7-ray emission of blazars indicates 
a double-peak structure in the overall spectral energy dis- 
tribution (SED) (Ulrich et al. 1997 and references therein; 
Fossati et al. 1998), suggesting two broad spectral compo- 
nents. The first component is generally interpreted as be- 
ing due to synchrotron emission, and the second is believed 
to be inverse- Compton emission of the same electron pop- 
ulation (e.g. Ulrich et al. 1997; Urry 1999). According to 
the different peak frequencies, BL Lac objects are divided 



into two subclasses, namely the low-energy-peaked BL Lac 
objects (LBL) which have synchrotron peaks in IR/optical, 
and high-energy-peaked BL Lac objects (HBL) which have 
synchrotron peaks in the UV/soft X-ray region. All of the 
above TeV 7-ray blazars are HBLs (the radio-selected BL 
Lac object 3C 66A has an arx = 0.74 < 0.80, Fossati et 
al. 1998) while almost all GeV 7-ray sources are LBLs and 
FSRQs. 

Although TeV 7-ray observations have already signifi- 
cantly affected our understanding of BL Lac objects and 
the unified scheme (e.g. Catanese & Weekes 1999), the 
sample of TeV blazars is still small. Much observational 
and instrumental work has been devoted to searching for 
new TeV blazars (Roberts et al. 1999; Urry 1999; Catanese 
& Weekes 1999 and references therein; Aharonian et al. 
2000). As Bai & Lee (2001) pointed out, according to the 
unified scheme of FRI radio galaxies and BL Lac objects, 
the SED of the jet-related emission of FRIs should also 
exhibit double-peak structure and some of them should be 
HBL-like. Though the relativistic boosting is week or even 
negative in FRIs, some nearby HBL-like FRIs may be loud 
enough to be detectable in the TeV energy range. If this 
turns out to be true, it will not only enlarge the sample of 
TeV sources but also provide a new approach for the study 
of radio galaxies. Moreover, TeV observations of nearby 
FRIs will provide a test of the unified scheme. 

Cen A (NGC 5128) is the prototype of FRI radio galax- 
ies, at a distance of 3.4 Mpc (z=0.0008) with two strong 
X-ray jets (e.g. Kraft et al. 2000). It is the only GeV 
7-ray source not belonging to the blazar class (Sreekumar 
et al. 1999). M87 (3C 274, Virgo A, NGC 4486, z=0.0043) 
is a well studied FRI radio galaxy with the brightest op- 
tical jet. It shows most of the characteristics of BL Lac 
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objects, with the jet oriented at 30° — 35° to our line of 
sight (Tsvetanov et al. 1998 and references therein). 

In this paper we present evidence showing that Cen A 
and M87 are HBL-like FRIs and that both of them are 
probably loud enough to be detectable at the TeV energy 
range during outburst. 

2. PROPERTIES OF 7-RAY-LOUD BL LAC OBJECTS 

The 7-ray-loud BL Lac objects (FSRQs similar to LBLs) 
have some common properties which imply physical sim- 
ilarities among different objects. It is believed that 7- 
ray emission from TeV BL Lac objects is dominated by 
the synchrotron self-Compton (SSC) process. Within the 
present uncertainties, the low frequency side of each of the 
two peaks of TeV BL Lacs can be described by the same 
spectral index (Tavecchio et al. 1998), consistent with the 
SSC model. Besides the property of double-peaked SEDs 
mentioned in the Introduction section, 7-ray-loud blazars 
have the following properties which should also be shared 
by HBL-like FRI radio galaxies, according to the unified 
scheme: 

1) The relation between the two peak frequencies 

In the context of the SSC model, the relation between 
the peak frequencies of the synchrotron component and 
the Compton component fc is Uc/fs oc 7peafc, where 7peafc 
is the characteristic electron energy (Urry 1999). For 
TeV BL Lac objects, the averaged upshifting factor of the 
Compton peak relative to the synchrotron peak is 10^^^, 
i.e., 



10 



,8±1 



(1) 



For example, Mrk 421 has a synchrotron emission peak at 
or slightly above lO^^Hz, and a Compton peak just below 
1 TeV (Urry 1999), yielding i^c/z^s <, 10^. The synchrotron 
peak of Mrk 501 occurs in the range 2 — lOOkeV (Catanese 
et al. 1997; Pian et al. 1998; Urry 1999; Catanese & Sam- 
bruna 2000), and the Compton component peaks at ~1 
TeV (Pctry et al. 2000; Sambruna 1999), with Vc/vs ^ lO'^ 
- 10^. PKS 2155 — 304 has a peak in its synchrotron com- 
ponent in the soft X-ray range (Edelson et al. 1995; Urry 
et al. 1997; Chiappetti et al. 1999) and a peak in its 
Compton component between lOGeV and 1 TeV (see Fig. 
5 of Chadwick et al. 1999), probably at lOOGeV (Kataoka 
et al. 2000; Chiappetti et al. 1999) with i^c/i's ~ 10^-10^. 

Since the flux sensitivity of TeV 7-ray detectors is not 
high, it is necessary for a candidate TeV source to have a 
peak in its Compton component in the TeV 7-ray energy 
range in order for it to be detected. BL Lac (2200-1-42, 
z=0.069) is one of the brightest LBLs. It has a Comp- 
ton radiation peak in the GeV energy range (Fossati et 
al. 1998) and is a strong GeV 7-ray source but it was 
not detected in the TeV energy range even during a large 
optical/GeV outburst in July/August 1997 (Bloom et al. 
1997; Bai et al. 1999; Aharonian et al. 2000). Thus, TeV 
candidates should be selected among HBLs or HBL-like 
objects, according to Eq. (1). 

2) The relation between the two peak fluxes 

The two radiation components of blazars are self-similar 
in shape. In the context of the SSC model, the luminosi- 
ties {Lg, Lc), peak luminosities {Ls{vs) and Lc('^c)), and 



peak flux densities (fai^s) and fc{vc)) of the two compo- 
nents are related by (Tavecchio et al. 1998; Stecker et al. 
1996) 

During the high state, the five TeV blazars all have a 
Compton luminosity Lc comparable to or slightly less than 
the synchrotron luminosity Lg, i.e., Lc/Lg ~1. Assuming 
this is also valid for candidate TeV sources, we obtain 



(3) 



According to the SSC model, the Compton component 

has the same spectral shape as the synchrotron compo- 
nent. Thus, with Eq. (3) and the spectral index of the 
synchrotron component, the TeV 7-ray flux of a candidate 
can be estimated as 
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where vi {nui < Vc) and V2 {1^2 > ^c) are the energy 
thresholds of a TeV 7-ray detector, and ai and a2 are 
spectral indices of the synchrotron component below and 
above Vg, respectively. 

3) The variability 

Perhaps the most striking property of TeV 7-ray sources 
is that the TeV flux is strongly variable. The simultane- 
ous multiwaveband observations for Mrk 421 and Mrk 501 
show that the amplitude of variations at the TeV energies 
is the largest (Retry et al. 2000; Catanese et al. 1997; 
Takahashi et al. 1996). Both of them exhibited low states 
in which the TeV flux is below the limit of detection. It 
is thus obvious that during the low state Lc/Lg <1, and 
sometimes L^Lg <<1. For this reason, Eqs. (3) and (4) 
are only valid during outburst. 

For Mrk 421 and Mrk 501, TeV flux varies with X- 
rays without a significant time lag (<24hr). This may be 
caused by changes in electron distribution, so it is likely 
that other TeV AGNs may also show similar behavior. 
Therefore, the TeV 7-ray outburst of an object can be 
predicted by investigating the X-ray variability of the ob- 
ject. 

4) Beaming effects 

For 7-ray-loud blazars, relativistic beaming is necessary 
not only to enable the 7-ray photons to escape from the 
source, but also to amplify the fiux and therefore make the 
source more easily detectable (e.g. Urry & Padovani 1994; 
GhiscUini et al. 1998). The observed {Fobs) and intrinsic 
{Fq) fluxes from a continuous relativistic jet are related by 



Fobs = S'+"Fo = 6PF0, 



(5) 



where S is the Doppler factor and a is the spectral index. 
The Doppler factor is defined as 5 = [r(l — Pcos9)]~^, 
where is the speed of emitting plasma in units of the 
speed of light, F is the Lorentz factor, defined as P = 
(1 — /3^)~^/^, and 6 is the viewing angle. It is notable that 
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beaming does not affect tfie reiations between two com- 
ponents discussed above (see Fig. 5 in Chiaberge et al. 
2000). 

Based on the TeV variabiiity, it was estimated that Mrk 

421 and Mrk 501 have 0^6° (Padovani 1999). If they 
were observed at 6* = 60°, according to the Eq. (5) and 
assuming typical values for F and p, i.e., F ^5 and p ^3, 
the observed flux would decrease by a factor of 1.3 x 10"'', 
which is not detectable even during the largest outburst. 
However, if they were at the distance of Cen A, z=0.0008, 
the observed flux would increase by a factor of about 
(§m^)^- Therefore, if Mrk 421 and Mrk 501 were moved 
to the position of Cen A and observed at a viewing angle 
of 60°, or moved to the position of M87 (z=0.0043) and 
observed at ^ = 35°, they would still be detectable at TeV 
energies during outburst. 

3. CEN A AND M87 AS TEV SOURCE CANDIDATES 

As mentioned in Section 1, BL Lac objects are divided 
into HBLs and LBLs according to the peak locations in 
the SEDs. Further studies show that for LBLs the X- 
rays are inverse-Compton emission, while X-rays in HBLs 
are an extension of synchrotron emission at lower energies 
(Padovani & Giommi 1996; Fossati et al. 1998). In the 
following subsections we will show that X-rays in Cen A 
and M87 consist of synchrotron emission and thus they are 
HBLs. 

3.1. Cen A 

Figure 1 shows spectral energy distributions of Cen A, 
M87 and Mrk 501 based on several observations in the 
literature. The X-ray spectral index of Cen A nucleus 
in 2-10 keV was found to be less than 1 (Turner et al. 
1997; Miyazaki et al. 1996), and hence \og{uf^) is a 
monotonously rising function of log i>. Between 10 and 
100 keV, the spectral index is 0.7 — 0.8 (Miyazaki et al. 
1996; Bond et al. 1996; Wheaton et al. 1996; Kinzer et 
al. 1995), also less than 1, with a monotonously rising 
SED of Cen A in this energy range. The SED of Cen A 
continuously rises below ^150 keV where the spectrum 
breaks, then the SED goes down continuously to the MeV 
and even GeV energy range (see Fig. 1). As in the case 
of blazars, this hump represents one of the two radiation 
components, and the radiation mechanism of this compo- 
nent is the same as that of X-rays. 

Steinle et al. (1998) and Kinzer et al. (1995) pointed 
out that the spectra of Cen A show interesting similar- 
ities to those of jet-aligned blazars (McNaron-Brown et 
al. 1995), and in particular to those of the well-studied 
quasar 3C 273 (.Johnson et al. 1995; Licliti (-t al. 1995); 
both show spectral breaks in the soft 7-ray regime, and 
both have intensity-independent power-law shapes below 
the break. For 3C 273 the peak in the soft 7-ray regime 
is a Compton peak, so Steinle et al. (1998) and Kinzer et 
al. (1995) regarded the peak of Cen A at ~150 keV as a 
Compton peak. 

However, based on the observations by ROSAT and 
ASCA, Turner et al. (1997) found that the spectrum of 
the jet of Cen A is actually consistent with the spectrum 
of the nuclear continuum, and that the spectrum of the 
jet is flatter than that expected as a result of an inverse- 
Compton scattering of radio photons but is consistent with 



that predicted by a simple synchrotron model. Recent 
Chandra observations also indicated that X-rays from the 
jet of Cen A comprise synchrotron emission (Birk & Lesch 
2000). Therefore, the peak at ~150 keV is not the Comp- 
ton peak but the synchrotron peak, as the case of Mrk 
501 (sec Fig. 1). The peak frequency of Cen A is variable, 
which is also similar to the behavior in TeV blazars. Obser- 
vations by CGRO OSSE in 1991 -1994 revealed peak fre- 
quencies of ~150 keV (Kinzer et al. 1995), while Miyazaki 
et al. (1996) reported a peak frequency of 180 keV, based 
on Welcome- 1 observations. 

The peak of the Compton component of Mrk 501 is at 
'^1 TeV, so that of Cen A is probably also at ~1 TeV. 
During the outburst state, the flux density at 100 keV is 
« 11 X 10"''' photon cm^^s^^keV^^ (sec Fig. 5 in Bond 
et al. 1996), and the synchrotron peak flux density is 
fi50kev = 2.1 X 10~^^erg cm~^s~^Hz~^, assuming a spec- 
tral index of 0.7. According to Eq. (4) with ai=0.7 and 
a2=1.3, the estimated TeV flux in 0.25 - 30 TeV for Cen 
A is 

F(0.25 - 30TeV) = 6.4 x lO-^gcm-^s"^ (6) 

which is about 10 times brighter than that of Mrk 501 in 
the 1997 outburst, suggesting that Cen A can be a strong 
TeV 7-ray source during outburst. 

3.2. M87 

The radio to X-ray spectral energy distribution of M87 
is very smooth like that of a typical BL Lac object (Tsve- 
tanov et al. 1998; see Fig. 1). If it is really a misaligned 
BL Lac object, the radio to X-ray hump should be the 
synchrotron component, and X-rays are also synchrotron 
emission, suggesting that M87 is an HBL. EINSTEIN and 
ROSAT observations show that the X-ray emission from 
the jet of M87 is probably the high frequency tail of the 
radio to optical synchrotron spectrum, although other ori- 
gins (thermal or inverse-Compton) for the X-rays cannot 
be ruled out completely (Biretta, Stern & Harris 1991; 
Neumann et al. 1997). The EUVE observations of M87 
are consistent with the spectral cutoff in the spectrum of 
the jet in M87 as suggested by Meisenheimer et al. (1996), 
which further supports the idea that the EUV and X-ray 
emission of the jet is synchrotron radiation (Berghofer et 
al. 2000). Reynolds et al. (1999) also identify M87 with an 
HBL, by combining Rossi X-ray Timing Explorer (RXTE) 
hard X-ray (4 - 15keV) observations with ROSAT data 
and inferring that the X-ray spectra of the M87 core and 
jet must be steep like those of HBLs. 

The synchrotron peak of M87 lies between 10'^''' — lO"'^ 
Hz (see Fig. 1). According to Eq.(l), the Compton com- 
ponent of M87 may peak at ~0.1 TeV. The flux of the 
core and jet of M87 in the EUVE DS band pass is 16.25 
/xJy (using a core/jet flux ratio of ~1.5) with a spectral 
index of the UV to X-ray power-law spectrum of a « 1.4 
(Berghofer et al. 2000; Harris et al. 1997). According to 
Eq. (4), the estimated TeV flux (0.25 - 30 TeV) for M87 
is 

i^(0.25 - 30TeV) = 1.1 x lO-^crgcm^^s^^ (7) 

which is comparable to that of Mrk 501 in 1995, suggesting 
that M87 may be a TeV 7-ray-loud FRI radio galaxy. 
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4. DISCUSSION AND CONCLUSIONS 

Cen A was observed with the CANGAROO 3.8m TeV 
7-ray telescope from March to April 1995, but was not 
detected (Rowell et al. 1999). There was no X-ray ob- 
servation for Cen A in 1995, but using the observations 
before 1995 we can infer that Cen A was in the low state 
during 1995. Between 1992-1995 Cen A was very faint and 
there was no trend towards outburst in 1995 (see Fig. 5 
in Bond et al. 1996). Thus, the TeV flux of Cen A was 
below the detection limit of the CANGAROO telescope. 
Grindlay et al. (1975) reported that Cen A was detected 
in the TeV energy range during 1972-1974. Since the de- 
tection has not been repeated (Israel 1998), it appears that 
no one has included Cen A in the list of known TeV 7-ray 
AGNs. However, from 1972 to 1976, Cen A underwent a 
large outburst in the X-ray energy range (see Fig. 4 in 
Turner et al. 1997 or Fig. 5 in Bond et al. 1996), and 
probably underwent an outburst in the TeV energy range 
at the same time, so the detection of Cen A by Grindlay 
et al. (1975) is considered to be authentic. 

At present, the CANGAROO 10m telescope is the only 
TeV 7-ray telescopes which can be used to observe the 
southern object Cen A. In the future, Cen A can be 
observed by CANGAROO III and the German-French- 
Italian experiment HESS in Namibia. The long-term X- 
ray light curve of Cen A shows that the period of outburst 
is 5 — 6 years (see Fig. 5 in Bond et al. 1996), which is 
typical for radio- loud AGNs (Smith and Nair 1995). The 
recent observation by the Chandra X-ray Observatory in 
1999 shows that Cen A is in the low state at present (Kraft 
et al. 2000). Thus, it is expected that Cen A may undergo 
an outburst in the near future. The long-term light curve 
also shows that Cen A did not undergo any large outburst 
since 1986, so the coming outburst may be a large one. In 
addition, according to Eq. (1) the Compton component of 
Cen A may also peak as high as ~100 TeV. It may thus 
be strong enough to be detectable even in the PeV energy 
range during outburst. 

The HST has recorded that the flux from the M87 opti- 
cal nucleus varies by a factor of ~2 on timescales of ~2.5 



months (Tsvetanov et al. 1998). A difference by a fac- 
tor of five in the flux detected by Spacelab 2 and Ginga 
indicated that the hard X-ray from M87 core also varied 
violently (Guainazzi & Molendi 1999), which suggests that 
the variability in M87 around peak frequency is similar to 
those in TeV blazars. EGRET has observed M87 but did 
not detect any GeV 7-rays (Sreekumar et al. 1996), which 
may indicates that M87 does not have a peak in its Comp- 
ton component at the GeV energy range, being consistent 
with our result. M87 was observed at TeV energies many 
years ago but was not detected. At that time the TeV 
detectors were not sensitive (the Crab Nebula was not de- 
tected either), and the observations only gave upper limits, 
F{> 0.4TeV) < 8.3xl0""photonscm-2s-i (Cawleyet al. 
1985) and F{> 0.21TeV) < 1.2 x lO^^Ophotonscm-^s-i 
(Weekesetal. 1972; see Fig. 1). M8 7 has not yet been ob- 
served by any modern TeV 7-ray telescopes. It can be ob- 
served using powerful TeV 7-ray telescopes, such as Whip- 
ple and HEGRA. 

In summary, the X-rays from Cen A and M87 are syn- 
chrotron emission rather than inverse-Compton emission, 
and hence these nearby FRI radio galaxies are HBL-like 
objects. In particular, Cen A is similar to Mrk 501 with 
its synchrotron component peaking at very high energies, 
which was once thought to be non-synchrotron radiation 
(Skibo et al. 1994; Kinzer et al. 1995; Steinle et al. 1998). 
According to the unified scheme of BL Lac objects and 
FRI radio galaxies and assuming that the properties of the 
known TeV BL Lac objects are common for all TeV AGNs, 
we predict that Cen A may have a peak in its Compton 
component power output at ~1 TeV, and that M87 may 
have a Compton emission peak at '~0.1 TeV, both having 
TeV 7-ray flux detectable by TeV 7-ray detectors available 
today. 
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Fig. 1. — Spectral energy distribution of Cen A (filled circles), M87 (asterisks) and Mrk 501 (open circles). For Cen A, the data from radio 
to optical bands are from Marconi et al. (2000), and X-ray and soft gamma-ray data obtained in 1991 are from Steinle et al. (1998). For 
M87, EUV data are from Berghofer et al. (2000), UV data are from Sparks et al. (1996), and other data are from Biretta et al. (1991). For 
Mrk 501, the SED is based on simultaneous observation of MJD=50556.32 (Petry et al. 2000). The upper arrow represents upper limit from 
Weekes et al. (1972), and the lower arrow represents upper limit from Cawley et al. (1985). 
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